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Abstract 

The one-dimensional extended isotropic XY model (s=l/2) in a transverse field with uniform 
long-range interactions among the z components of the spin is considered. The model is exactly 
solved by introducing the gaussian and Jordan- Wigner transformations, which map it in a non- 
interacting fermion system. The partition function can be determined in closed form at arbitrary 
temperature and for arbitrary multiplicity of the multiple spin interaction. From this result all 
relevant thermodynamic functions are obtained and, due to the long-range interactions, the model 
can present classical and quantum transitions of first- and second-order. The study of its critical 
behavior is restricted for the quantum transitions, which are induced by the transverse field at 
T = 0. The phase diagram is explicitly obtained for multiplicities p = 2, 3, 4 and oo, as a function of 
the interaction parameters, and, in these cases, the critical behavior of the model is studied in detail. 
Explicit results are also presented for the induced magnetization and isothermal susceptibility Xt-> 
and a detailed analysis is also carried out for the static longitudinal {S^Sf) and transversal {SJSf) 
correlation functions. The different phases presented by the model can be characterized by the 
spatial decay of the these correlations, and from these results some of these can be classified as 
quantum spin liquid phases. The static critical exponents and the dynamic one, z, have also been 
determined, and it is shown that, besides inducing first order phase transition, the long-range 
interaction also changes the universality class the model. 

PACS numbers: 05.70.Fh\sep 05.70.Jk\scp 75.10.Jm\sep 75.10.Pq 

Keywords: One-dimensional extended XY-model; long-range interaction; quantum phase transitions; static 
properties; quantum spin liquids. 
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I. INTRODUCTION 



The study of the quantum phase transitions(QPTs) has been the object of many theoret- 
ical and experimental investigations [1], since they play an essential role in understanding the 
low temperature properties of the materials ^ . In particular, for low dimension magnetic 
materials, the study of quantum spin chains has provided much insight in this direction^. 
Motivated by this, in the last years, others tools have been used to investigate QPTs in spin 
chains, such as, and quantum fidelity [71], where 



the last two concepts were unified in the approach of the geometric tensors[8]. All these tools 
have only one purpose, that is, to point out the existence of the quantum critical behavior. 
On the other hand, the QPTs are also directly related to quantum entanglement, which 
has been object of study in quantum computation [9]. Therefore, the study of spin systems 
is of great importance for understanding the behavior of the materials in low temperature 
regime. 

Among the quantum spin models, the one-dimensional XY model introduced by Lieb et 
a/, [h]], with its generalized classes, have attracted much interest in the last decades. In 



particular, the models with multip 



long-range interactions(see Refs. [13 



e spin (see Refs. 



11 



12], and references therein) and 



141 ]. and references therein) have been object of inten- 
sive investigations. As pointed out by Derzhko [l^] and in the references therein, quantum 
spin systems with multiple spin interactions work as effective spin models for the standard 
Hubbard model under certain conditions. On the other hand, models with long-range inter- 
actions are important to understand the process of quantum information in spin chains as 
well as for the study of the classical and/or quantum crossover, as it was explained by de 
Lima and Gongalves.[3]. Another importance of the long-range interaction is that it can 
induce first-order QPTs, which play an important role in the quantum critical behaviour[l3]. 

Therefore, in this paper, we will study the effect of the long-range interaction on the 
quantum critical behaviour of the extended one-dimensional XY-model with multiple spin 



interactions 



181 ] . The solution of the model can be obtained exactly, at arbitrary temper- 



atures in a one-dimensional and the main purpose of the paper is to study the quantum 
critical behaviour of the model. 

We analyze these two kind of interactions, where we will show the role of each one, that 
is, we will show that the long-range interactions induce first and second order QPTs and 
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that their presence change the universality class of the model, as we will also verify the 
scaling relations proposed by Continentino and Ferreira[l^] for first-order QPTs. While for 
the multisite interaction among spins, we will show the presence of many different kinds of 
quantum spin liquid phases, which after some time returned to be a relevant topic in the 
present research [20]. 

In Sec. II we introduce the model and obtain its exact solution by means of Jordan- 
Wigner fermionization and the integral Gaussian transformation, and present its basic 
results, such as, the spectrum of the energy, magnetization at arbitrary temperatures. The 
quantum phase diagrams, as function of the interaction parameters, are presented and dis- 
cussed in Sec. Ill, and in Sec. IV we study the scaling behaviour of the longitudinal and 
transversal correlations functions on the different phases presented by the model. In Sec. V 
we evaluate the critical exponents, and discuss the change in the universality class due the 
presence of the long-range interactions. Finally, in Sec. VI we summarize the main results. 



II. THE MODEL AND BASIC RESULTS 



We consider the one-dimensional isotropic extended XY model (s = 1/2) with uniform 
long-range interactions among the z components of the spins, whose Hamiltonian is explic- 
itly given by 

Y j S S X 

N 



N j N N N 

H = - Ji ^2(Sj Sj +1 + SjSj +1 ) - — ^2 S j s i ~ S j~ 



5=1 
N V 



j=l 1=1 



3=1 



T \ ^ \ ^ ( qX qX , qy qy \qZ qZ qZ 

j=\ k=2 



(1) 



where the parameters J\, J2 are the exchange coupling between nearest neighbors, I the 
uniform long-range interaction among the z components, p is the multiplicity of the multiple 
spin interaction, and iV is the number of sites of the lattice. 
Introducting the Jordan- Wigner transformation 



oxi> I m c ] c l 
1=1 



c] ; S j - Cj 



°j ~ c j c j 2 ' 



oxp ( —in cjci 
1=1 



(2) 
(3) 
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where c, and c] are 



ermion operators, the Hamiltonian can be written in the well known 



decomposed form( 211] also references therein) 

H = H + P + + H~P- 



where 



N N 



3=1 3=1 3=1 1=1 

- f E E (t) ( c ^« + c 5« c ') + t (* - 

3=1 fc=2 V 7 V 7 

± y ( c !v c i + C W) ± y ^ ( ~2~ ) ( c lv- P +« c ;-P+fc + cL P +fc c w-P+«) ± 



(4) 



(5) 



1=1 k=2 
p-l>2 



±fE E (t 

k>2 



and 



with P given by 



( N \ 
P = exp j «7T c|q J . 



(6) 



(7) 



As it is also well known, the Hamiltonian 7i ± can be diagonalized by imposing antiperi- 
odic (for 7i + ) and periodic (for H~) boundary conditions, and, in the thermodynamic limit, 



the static properties can be described by 7i~ 



21 



22 



231 ] . Therefore, since we are interested 



in the determination of the static properties, we will identify the Hamiltonian of the system 
with 7Y~. By taking into account that the long-range interaction term commutes with the 
Hamiltonian, the partition can be written in the form 



Zj^f = exp 



2 V 2 

JV 



Tr <^exp 

N p 



BJ N 

~Y 2 ( c ^'+i + 4+i c i) + 



3=1 

K-l 



+ P (h - I) £ c)c, + § £ £ (± ) + cj + . 

3=1 3=1 K =2 V 7 



xexp 



31 N N 

j^J2J2 c h c lci 

3=1 1=1 



(8) 
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where (3 = 1/ksT and T is the temperature. Introducting the Gaussian transformation 24 ] 



exp{a 



exp ( — — + V2ax ) dx, 



the partition function can be written in the integral representation as 



(9) 



Zx = exp 



^(s_I 

2 V 2 



^ f exp ( — \ lr \ exp 




— oo 

N -J N p , s k-1 

+ (5 - 7 + 2 Vre) X) =>i + y E E (t) 



2 



-jCj +K + Cj +K Cj 



( c h+ i + c ]+i c i) + 

(10) 



3=1 



dx, 



where x = x/VN, J\ = (3J\, J 2 = (3J 2 , h = (3h and I = (31. 
Introducing the canonical transformation 

N N 



—L exp (ijq) c q and c] = -== ^ exp (-ijq) cj, 

* (7=1 * 0=1 



(11) 



9=1 9= 

with g = 27m/ N, n = 1, 2, AT, Eq. fllOp can be written in the form 

Nx 2 * 



where 

cm 

with 



N 




27T 



exp 



Zu = C(/3) / exp 



2 V 2 



(12) 



C (a;) = Tr < exp 



E 



6 q [X) C' q C q 



€ q (x) = (3 



Ji cos(g) + J 2 ^2 ( ~~ o ) cos(Kg) + h — I + 



K=2 




X 



(13) 



(14) 



In the thermodynamic limit, N — » oo, we use the Laplace's method 25] to evaluate the 
partition function, and Zj^ can be written in the form 



exp 


N 
2 




+ Ng(x ) 




g"(x ) 


1/2 



where 



with g(x~o) satisfying the conditions 

g'(x ) = 0, and g"(x ) < 0, 



(15) 



(16) 



(17) 
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and xq explicitly given by 



Ti N 



Xq 



N 



^ 1 + exp(-e g (x ))' 



(18) 



Therefore, from Eq. (1151) we can write the Helmholtz free energy as 

F N = y (h - - Nk B Tg(x ) + *|L In | /(x ) | . (19) 
Taking into account that the magnetization per site M z can be written in the form 



M z 



1 



1 



2 iV^l + exp[-e 9 (x)] 2 : 



j i q 

by using Eq. (TTB"j) , we can express xq in terms of M z in the form 

1 



x 



21 M 



(20) 



and from this result it follows that the functional of the Helmholtz free energy per site is 
given by 

/ = J + IM Z (1 + M z ) + ^ f In (1 + exp [e q (M z )}) dq, (21) 
2 ^ Jo 

with 

_ _ p / lV -1 

e q (M z ) = -J lC os(q) - J 2 Y] (--) cos(Kq) - h-2IM z . (22) 



K = 2 



From Eq. fl2T]) . we can determine the equation of state imposing the stability conditions 

5/ 



<9M 22 



= 0, 
> 0, 



which leads to the result 



2tt 



/ tanh 




/o 


2 



(23) 
(24) 

(25) 



III. QUANTUM CRITICAL BEHAVIOR 

In the limit T — > 0, the functional of the Helmholtz free energy per site[Eq. fl2T]) ] is given 



by 



/ = £+IM*(l + M*)-- / e g (M z )dq, 
2 71 J 



(26) 



e,(M z )<0 



where 

K-l 



e q (M z ) = Jicos(g) + J 2 ^ ( -- J cos(«;g) + h + 2/M 2 , (27) 

and the equation of state, given by Eq. (|25|) . can be written in the form 

M 2 = 7 1 sign[e q {M z )}dq. (28) 

The quantum phase diagram of the model, for second order phase transitions and arbi- 
trary p, can be obtained from the previous expression by imposing the divergence of the 
isothermal susceptibility, x¥ = dM z /dh — > oo, and for first-order phase transitions, by 
using the equation of state, Eq. ( 1251) . and by imposing the condition 



/(M 2 ) = /(MJJ, (29) 

where M? t and ML are the values of the induced magnetization at the transition. 
Following Titvinidze and Japaridze ll| , by introducing the unitary transformation 

gx,V _ s *,v . gz ,S ; . (30) 

and the time reversal one 

qX,Z qX,Z . qV _ qV I o 1 \ 

ij • ' ij • ' i — i ' v / 

we can show that the Hamiltonian is invariant under the transformation J\ — > — Ji and 
J 2 — > — J 2 - This means that only the signs of 7 and J 2 are relevant in determining the critical 
behaviour of the system, and the appearance of multiple phases are result of the competition 
between these interactions which induce frustration in the system. In particular, they can 
induce the so-called quantum spin liquid phasesf^], as in the case of the model without long- 



range interaction [HI]. Therefore, without loss of generality, we will consider in all results 
presented J\ > only. 

From Eqs. f!26H28p . we can show that, for arbitrary p, the system presents second order 
quantum transitions for / < 0, and first-order quantum transitions / > 0, as in the previously 
studied XY-models with similar long-range interactions [l^]. 

Although the Eqs. ( 12611281) allow us to obtain the phase diagrams for arbitrary p, we 
will only consider the cases p = 2,3, 4 and oo, since they present the main features and, in 
these cases, analytical expressions can be obtained for the critical surfaces and critical lines 
associated to the second order quantum transitions. 
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For the case p = 2, which has also been studied by Titvinidze and Japaradze [1 1| and 
Krokhmalski et al. |12J for the model without long-range interactions, we show in Figs. Q] 
and [2] the critical field h/J\ as a function of I / J\ in the regions J2/J1 > and J2/J1 < 0, 
respectively, for different values J2, which are projections of the global phase diagram. It is 
worth mentioning that in this case the results for J%j J\ < can be obtained from the results 
for J2/J1 > by introducing the transformations h/J\ — > —h/J\ and J2/J1 — » —J2/J1, a s 
can be verified in the results shown in Figs. [1] and [2j 

In Figs. |3] and H] we present the magnetization and the isothermal susceptibility, respec- 
tively, and from their behavior we can conclude that for I / J\ < the system undergoes 
second order transitions, and first-order transitions for I / J\ > 0. 

As can it be seen in Figs. [T] and [2J the number of transitions of first and second or- 
der depends on J2, and it can also be shown that these transitions correspond to three 
phases for < J2/J1 < 1/2, and to four phases for Jij J\ > 1/2. Following Titvinidze and 
JaparandzeQ, we can classify the intermediate phases, whjch are limited hy the two satn- 
rate ferromagnetic phases, as quantum spin liquid phases[20]. As will be shown later, these 
spin liquid phases will be characterized by the spatial decay of the transversal static corre- 
lation function (SjSf) and by the modification of the oscillatory modulation of longitudinal 
static correlation function {S^Sf). 

The global phase diagram for p = 2 is shown in Figs. [5] and [6], for J2/J1 > and 
J2/J1 < 0, respectively. For / < 0, the critical surfaces can be obtained explicitly and, for 
J2/J1 > 0, where there are four critical surfaces, which are given by the equations 

k-i-b^'^T^' (32) 

_^_|_ 1 .0,lor0<£4 (33, 
h J\ I h J 2 1 , . 

-2T 1 -U2 + T 1 -T 1 = ^ ioT J1-2' (34) 

j 2 21M? h j 2 1 . . 

^- — -T^^^^T^r (35) 

where 

arccosf^ 1 — 1) 1 

Ml = ^ - - (36) 

7T 2 
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and for J 2 j J\ < there are four critical surfaces given by 



Jo I h Jo 

J 2 I ^ + ! _ for 1 < J<2 < 

2,J\ J\ J\ 2 J\ 

J 2 J\ I h J 2 1 



J 2 2IM, 



h 



Jo 



(37) 

(38) 
(39) 
(40) 



where 



Mi 



arccos 



■1l 



+ 1) 1 



71 2 

The critical lines, shown in Figs. (HandEJ can be obtained from Eqs 
and J2/J1 = 2.0, and from Eqs. EMU for J2/J1 = -0.2 and J 2 /Ji = -2.0. 

It should be noted that for J 2 j J\ = 1/2 the critical surfaces meet at a bicritical line[l3j 
given by 

h 3 T 

(42) 



(41) 

or J2/J1 = 0.2 

n 



A 



3 1 



For / > 0, where the phase transitions are of first-order, critical surfaces are obtained 
numerically from the solution of the system 



0. 



M, 



fMZ^U + 1 = 0, 

7T I ' 



(43) 



where f 2t an d fit are given by 



fit 



fit 



arccos 



arccos 



f(Mit) - f{Ml) = 0, 



Ji + VJ! + gj + 4 J 2 (g + 2/M/, 
2J, 



2J 2 



(44) 



(45) 



Mf t , Mj t are the values of the magnetization at the transition, with i, j = 1, 2, 3, 4 for 
< J2M < 1/2, and with z, j = 1, 2, 3, 4, 5, 6 for J 2 /J 1 > 1/2, and /(M|) is given by 
Eq. f[2"6T) . As in the case of second order transitions, the critical lines shown in the Fig. [TJ 
are determined from the previous systems by considering J 2 j J\ = 0.2 and J 2 j J\ = 2. 
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For J2/J1 > 1/2, there are four critical surfaces and three of them, for J2/J1 = 1/2, meet 
at a critical line which is determined by the following system of equations 



M- 



it 



M. 



21 



7T ~ 2 

fit- fit I 1 
7T ~ 2 



Ml - 



+ 



: 0, 

0, 
0. 



(46) 



}{M z lt ) - f(M§ t ) 
{ /(MfJ - /(M| t ) 

where / (M" t ) is obtained from the Eq. fl26l) . and with i = 1, 2, 3 and M| t = —1/2. This is 
a triple line which meets the bicritical line at / = 0. 

A second triple line can be determined by imposing M| t = — 1/2 and M% t = 1/2 in the 
system 

( ..2 

= 0, 



ML - ^ 

Zt tt 



/(M| t ) - /(Mf t ) = 0, (47) 
/(M| t ) - f(Ml t ) = 0, 

where /(M| t ), and v?2t are given by Eqs. (|2"S"|) and (145} . which begins at the point J2/J1 = 0, 
h/J\ = and J/Ji = 4/7r, which has been obtained Gongalves et al. [131 ]. 

For the special case p — » 00, we can find the critical surfaces by using the same procedure 
used in the case p = 2. In this case, due to many intersections of the critical surfaces, the 
global phase diagram becomes too complicated, as we will show below. Therefore, we will 
present some projections of the global diagram which contain the main characteristics of 
this diagram and are shown in Figs. [TlfTUl 

In this case, the fermion excitation spectrum, obtained from Eq. (1271) . is given by 

2 cos (2q) + cos (q) 



e q (M z ) = J x cos(g) - J 2 



+ h + 21 M z 



(48) 



5 + 4 cos (q) 

and from this result we can determine the equations of the critical surfaces for I j J\ < 
and J2/J1 > 0, which are 



J_2 
Jl 



I 

Jl 



h 

Ti 



3t/i J\ J\ 




1 = 0, for ^ > 0, 
Ji 

= 0, for < ^ < -, 
~ Ji~ 23' 



I_A = for^>- 
J\ J\ J\ 23 



(49) 
(50) 
(51) 
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J 2 21 M? h „ J 2 27 

1 " 0, for > -, (52) 



3</i t/i t/i J\ 23 

where 

arccos f jH 2 ~^0 1 
M{ = V^z^ll _ I. (53) 

7T 2 

Identically, we can show that for I / J\ < and J2/J1 < 0, the critical surfaces are 

*-l-.| + 1 = 0, fa -l<£<0, ,54, 

A + l_|_ 1=0 , for -3<^<0, (55) 



5.7,-9.7, A / /, 0, for — < — !-, (56) 



where 



4J1 y Ji v J i / J i J i J i ~ 11 

J 2 2/M| /i J 2 1 

j 1 "^r _ x +1=0 ' for ~ 3 - x-"H' (57) 

^ + 4-| + l = 0, for^<-3, (58) 
J± Ji J\ J\ 

3^-^r-^- 1 = ' for ^^- 3 ' (59) 

Ml = \ 4Jl 4J2; - ± (60) 

arccos ( I'V'+'y 1 ) , 

1 ; _ 1 (61) 

7T Z 

In this case there are three bicritical lines which are given by 

at J 2 jJ x = 27/23, (62) 

at J 2 /J x = -1/11, (63) 
h I 

— = -2 + — , at J2/J1 = -3. (64) 
J\ J\ 

For I / J\ > 0, as for the case p — 2, the first order transition surfaces can be determined 
numerically and the triple lines are determined by following the same procedure adopted for 
p = 2. 

In Fig. [7|the phase diagram is shown for the positive region J 2 / J\ > 0, where there are 
four critical surfaces which are given by the Eqs (I49H53I) . These critical surfaces meet at a 
bicritical line, which is given by Eq. 



h 


14 










Ti~ 


~23 


+ X 


h 


10 




T x ~ 


IT ~ 
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The phase diagrams for the negative region J2/J1 < and for different values of J 2 / J\ 
are shown in FigsJHlfTOl In this case there are six critical surfaces, given by Eqs I54"ll59l which 
meet at bicritical lines given by the Eqs.fJMKID- As we can see m FigsJSlTTUl the system 
presents in this case identical critical behavior to the one obtained for the case p = 2, as 
far as the critical behavior is concerned. However, it is worth mentioning the appearance of 
quadruple point at h/J\ = 0, J 2 /Ji = —1.7370..., and J/Ji = 1.7798...,which is shown in 
Fig. [9] and is not present in the case p = 2. The behavior of the functional of the Helmholtz 
free energy at this point is presented in Fig. [TTJ 

We have also analysed the cases p = 3 and p = 4, for I = 0, where the model presents 
new quantum spin liquid phases. Although in these cases there are no first-order transitions, 
we have restricted these analyses to the model without the long-range interaction, since the 
main purpose was to study appearance of new quantum spin liquid phases, which is mainly 
controlled by the multiple short-range interaction. 

In Fig. [12] we show the phase diagram for the case p = 3, where the critical lines are 
given by 

h 1 /13J 2 -12J x \i 19J 2 + 36Jx J 2 12 J 2 

— = , for — > — or — < —4, 67 

Ji 108 V Ji J IO8J1 ' Ji ~ 13 Ji ~ ' v ; 

h 1 /13J 2 -12Ji\^ 19J 2 + 36Ji J 2 12 J 2 

— = ,for — > — or — < 0. 08 

Ji 108 V Ji J IO8J1 Ji ~ 13 Ji v ; 

As it can be seen, there are six quantum spin liquid phases which, as we will show later, 
can be classified in three different classes as far as the critical behavior is concerned. 

The phase diagram for p = 4 is shown in the Fig. [13], where the critical lines are given by 

h 3J 2 . . 

rrsfr 1 ' (69) 

* " 2 +l, (70) 



Jl 8Jl 

h J 2 fRiV 3J 2 fRA 2 /5J 2 l\fR 1 \ 51J 2 1 



Ji 16 Ji V 2 / 32 Ji V 2 J V.I6J1 2 J V 2 J 256 Ji 4 
h J 2 /i? 2 \ 4 3J 2 [R 2 \ 2 fl 5J 2 \/i? 2 \ 51J 2 1 



(71) 



J L 16 Ji V 2 ) 327, \ 2 J + I 2 16jJ 1 2 J 256 Ji 4' ( ^ 72 ' ) 
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h 



Ji 8J1 



+ 



5J2 
8J1 



3; 



(■ 



5Uo 1 



V 256Ji 4 



where R\ = Ki+K 3 , with 



i? 2 = 5 3 + S 1 3 , with 



(73) 



K = (8J1 - 5J 2 )/J 2 + ^(64^ - 80 J X J 2 + 24 J|)/J 2 2 , 

J' = 166/(6 2 + 106 + 1) and b = [(3 + >/5)/2l 3 for J' < J 2 / ' J x < 4/3, (74) 



S = (5J 2 - 8J X )/J 2 + y (24Jf - 80 J X J 2 + 6472)771, for J 2 /J x < or J 2 /Ji > 2, (75) 
and i? 3i = cos {6 + 0j) /2, with 



9 = arccos [(8J1 - 5 J 2 )/ J 2 ] /3 and & = 2vr (i - 1) /3, i = 1, 2, 3 for 4/3 < J 2 /J x < 2. (76) 

In this for p = 3, the system presents six quantum spin liquid phases which can 

also be classified in three classes as far as the critical behavior is concerned, as well as there 
exist only second order phase transitions. 



IV. STATIC SPIN CORRELATIONS 



The static correlation function (SjSj +r ), in the thermodynamic limit, can be given by 



22 



21 



Tr \exp(-[3H-)S*SL r ] 

1 qz qz \ 1 f\ ^ ) j j+rj (77\ 

Tr[exp(-PH-)} ' [77) 
where 7i~ is the Hamiltonian given in the Eq. (j3J). 

After the introduction of the Fourier transform[Eq. ffTTl) ]. Hr can be written in the form 

^ = EM^)ch + y, (78) 
q 

where e q (M z ) is given by Eq. fl22|) . 

By introducing the fermion operators, we can write 

s i = jfY, e *pW«-fltt c *-l> ( 79 ) 

qq' 



14 



and from Eq. (!79|) . by using the Wick's theorem[26|, the static correlation can be written as 



(SjSj +r ) 



11 



with 



\ n i) 



1 



(80) 



(81) 



1 + exp [e q (M z )] 

From Eq. (|80|) . we can obtain the static correlation function which, after some straight 
forward calculations, can be written as 



(SjSj +r ) — 



-iVtanhf h. 



cos(gr) tanh ( ^ 



-i 2 



6, 



vo 



(82) 



and, in the thermodynamic limit, in the form 

2 

1 ? . /e 



-J tanh^Jrfg 



^/co s( . )tanh |H, 



<5r0 
4 ' 



(83) 



As it is well known 271 ]. at T = 0, the direct longitudinal correlation function of the short- 
range XY- model behaves asymptotically as p zz (r) = (S*S* +r )-(M z y ~ /(r)r" 2 , where /(r) 
is an oscillatory function. A similar behavior can be found for the above expression [Eql83] 
for the so-called quantum spin liquid phases, where the oscillatory factor f{r) undergoes 
changes for different phases, as shown by Titvinidze and Japaradze[ljJ], for the case p = 2 
without long-range interaction. 

In the presence of the long-range interaction, for p = 2 and J/Ji < 0, we have at most 
two multiple quantum spin liquid phases in the limit T — > 0. For the first one phase the 
parameters satisfy the conditions 



Jo I h Jo 
-1 + — + —< — < 1 + — 
2,J\ Ji Ji 2,Ji 



or 



(84) 



Ti 



-l + 



A 

2Ji 



Ji ~ 2' 

A^l< A<i + A 



(85) 



Ji' 
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where M( is given by Eq. (136|) . and from Eq. (I83I) we obtain 

sin 2 [gir] 



where 



P zz (r) 



qi = arccos 



7r 2 r 2 



Ji-y {M z 

2 J 2 



with 



y {M z ) = yJf + 2J| + 4J 2 (h + 21 M z ). 
For the second quantum spin liquid phase the parameters satisfy the conditions 



(86) 



(87) 



(88) 



J 1> 1 
Ji - 2' 

J 2 J\ I h J 2 2IM Z 



(89) 



2Jx 4J 2 Ji ~ Ji ~ 2Ji Ji ' 

and, in this case, the asymptotic behavior of the direct correlation p zz (r) is given by 



p zz (r) ~ - 
where gi is given by Eq. ( !87l) and g 2 by 

g 2 = arccos 



{sinfgir] — sin[g 2 r]}' 
ir 2 r 2 



Ji + y {M z 

2J 2 



(90) 



(91) 



Therefore, p zz (r) shows a similar behavior to the case where the model does not contain 
long-range interaction 11 ] . 



The correlation length which diverges at the critical point 28J , when we have a single 



phase, is associated to the period of t 
extension of its scaling form given by 



re oscillation of p zz (r) and is defined by an analytical 



22 1 



p zz {n) 



F(in/0 



(92) 



When the system presents two phases between the saturated ferromagnetic phases, be- 
sides the correlation lengths £i there is an aditional one £ 2 , associated to the adjacent 
transitions and, in this case, we can define a further extension of the scaling form as 
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and define lenght of the system £ in the form 



£ = & + (94) 

Therefore, for the case where we have a single phase, by using Eq. (|92|) . we can find from 
Eq. (|S5j) 



1 



2 arccos 



J\ — y {M z ) 



2 Jo 



(95) 
•1 + 



which diverges at the critical points given by h/J± = 1 + J 2 /2J\ — //Ji and /1/J1 = 
J 2 /2Ji + J/Jx. 

On the other hand, for the case where we have two phases between the saturated ferro- 
magnetic phases, by using Eqs. fl93|) . (JMj) and, by using Eq. fl90l) we can find 



1 



arccos 


" Ji+y(M z ) 
2J 2 


+ arccos 


" Ji-y(M*) 
2J 2 




2 arccos 


Ji-y{M*) 
2J 2 


n 


arccos 


2j 2 




— arccos 


Ji-V(M*) 

2j 2 


]} 



(96) 



which, as expected, diverges at the critical points h/J\ = — 1 + J2/2J1 — 2/Mj 3 /J 1 , h/J\ = 
1 + J2/2J1 — //Ji and /i/Ji = — J 2 /2Ji — Ji/4J 2 + //Ji, where AF is an intermidiate 
magnetization. 

Although the results presented are for p = 2, it can be shown that the correlation lenght 
always diverges at all second order critical points irrespective of the value of p, and for 
multiple transitions an aditional correlation length has to be introduced for each new phase 
presented by the system between the saturated ferromagnetic phases. Consequently, the 
scaling relation, given in Eq. (1931) . and teh correlation length of the system have to be rede- 
fined accordingly. 

The transversal correlation function (SjSj +r ), given by the Toeplitz determinant jlo| 

( (A.B2) (A^s) {A 1 B A ) ... (^1%*.+!)^ 
(A 2 B 2 ) (A 2 B 3 ) (A 2 B 4 ) ... (A 2 B j+r+1 ) 
(A,B 2 ) (A 3 B 3 ) (A 3 B 4 ) ... (A 3 B j+r+1 ) , (97) 



1 

4 



{(AjBj^) (AjB l+3 ) (AjB l+4 ) ... (AjB j+r+1 ) J 
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where 

(AjBt) = ~j cos[g(j - 0] tanh (^^-) dq, (98) 
o 

can be evaluated numerically. This correlation, for the usual short-range XY-model and at 
T = 0, behaves asymptotically as (SjSj +r ) ~ 27]. As in the case of the longitudinal 
correlation, the asymptotic behavior of the transversal correlation function (SjSj +r ) also 
presents changes in its power law decay due to the presence of the multiple spin interaction. 
This result has been shown by Titvinidze and Japaradze H[ , in the case p = 2 without 
long-range interaction, which classified these phases as spin liquid phases. 

For the model with long-range interaction, the transversal static correlation, at T = 0, 
was evaluated numerically from Eq. fl97j) by considering the maximum value of r equal to 
250, and the power decay determined by considering the scaling form (S^Sf) ~ /(r)r~^. The 
results for p = 2, are presented in Figs. [T41 and [151 where the scaling forms are shown in the 
insets. The transversal correlation function for the quantum spin liquid I phase, presented 
in Fig. [TH the scaling behavior is given by (SjSj +r ) ~ r~^, while for the quantum spin 
liquid II phase, presented in Fig. [15], it behaves as (SjSj +r ) ~ /(r)r _1 , where f(r) is an 
oscillatory factor. These behaviors are identical to the ones obtained for the model without 



long-range interaction 



111 ] .We would also like to point out that identical results are obtained 



for the case p —>■ oo. This general result, which depends on the number of phases only, gives 
support to the classification of the intermediate phases as spin liquid phases. 

The transversal static correlation was also calculated for p = 3, where there is a new 
quantum spin liquid phase, which can be identified in the phase diagram shown in Fig. [121 
In this new phase, denominated quantum spin liquid III phase, its scaling behavior is given 

3 ■ 

by (SjSj +r ) ~ /(r)r~2, where the oscillatory behavior /(r) is shown in Fig. [161 Finally, for 
p = 4, in the new phase denominated quantum spin liquid IV phase and shown in the phase 
diagram presented in Fig. [T3], the scaling behavior of the transversal correlation function is 
given by (SjSj +r ) ~ f(r)r~ 2 , which is presented in Fig. [T71 

It should be noted that, differently from the transversal correlation function, the spatial 
decay of the longitudinal correlation function does not depend on p. 
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V. CRITICAL EXPONENTS 



The critical exponents, at T = 0, associated with the magnetization, isothermal suscep- 
tibility, correlation length and the dynamic critical exponent z, for p — 2, can be evaluated 
analytically since the quantities of interest are known in closed form. Since these exponents 
are associated to second order transitions, our analysis will be initially restricted to / < 0. 

Therefore, following [29J, we define the order parameter given by 



\I = M z t - M z , 



(99) 



where M t 2 is the magnetization at the transition. This order parameter goes to zero at 
the second order transitions, and it is different from the one proposed by Titvinidze and 



Japaridze 



Then, from Eqs. ([28]), (JHTJ) and (JSHJ), we get 



M z = — arccos 

7T 



■h-y {M z 
2 J 2 



(100) 



and by expanding Eq. f llOOp up to second-order in M z , we obtain 

_2 



^(M^) 2 ^ Jl{hc - h \ 
2 V / Ji + 2J 2 



and 



M z 



Ji {h c - h) 
21 



for — 
Jl 



for — < 0. 
Ji 



0. 



(101) 



(102) 



From the Eqs. ( 11011) and ( 1102[) . and the scaling form M z ~| h c — h I 13 , we can conclude 
that the critical exponent (3 is given by /3 = |, when I / ' J x = and p = 1 for I/J 1 < 0, 
respectively, showing that the universality class has changed with the presence of the long- 
range interaction. 

The isothermal susceptibility can be obtained from the Eq. (11001) . and is given by 

Ji 



Xt 



nJ[y(M z )f 



1 - 



2J 2 



(103) 
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In the critical region, we have Xt H h c — h | 7 , and from the previous expression we find 
that the critical exponent 7 is equal to \, when ljJ\ = 0, and it is zero for I / J\ < 0. 
Since at T = 0, 7 is identical to a, we can show that in both cases, namely, with or without 
long-range interaction the exponents P, a and 7 satisfy the Rushbrook scaling relation 28]. 
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The critical exponent v associated with the correlation length £ ~| h c — h \ v , can be 
obtained from Eqs. fl95|) and (196]) . From these expressions we can immediately show that v 
is equal to |, when I/Ji = and it is 1 for lj J\ < 0. 



From the above results, by assuming the quantum hyperscaling relation 30] 



2-a = v(d+z), (104) 

where d is the dimension of the system and z the dynamic critical exponent, we find z = 1 
for / < 0, and z = 2 for / = 0. Therefore, we can conclude that the system also presents a 
non-universal critical dynamical behavior. 

For / > 0, following Continentino and Ferreira[lj|, we introduce a critical exponent 
associated to the free energy for the quantum first-order phase transition. Assuming the 
scaling form f(h ) = f{h t ) ± E ± \ h t — h \ 2 ~ a close to the field of transition ht, since the 
free energy given in the Eq. ( 1261) can be written as 



n , h + , „ , /if , 1 , h — ht , 

/(— ) = /(-^) - - | |, for h > ht, (106) 



./i t/^ 2 J\ 



we obtain a = 1, which gives support to Continentino and Ferreira conjecture 19] 



VI. CONCLUSIONS 

In this work we have considered the one- dimensional isotropic XY model with multiple 
spin interactions and uniform long-range interactions among the z components of the spin, 
in a transverse magnetic field. The solution of the model was obtained exactly for arbitrary 
p, which characterizes the range of multiple spin interaction, and at arbitrary temperature. 
Explicit equations have been obtained for the functional of the Helmholtz free energy from 
which the equation of state can be determined numerically. The quantum critical behavior 
was studied in details for p = 2, 3, 4 and oo, and the multiple phases presented by the model 
have been characterized by the asymptotic behavior of the static correlation functions 
(SjSf) and (SjSf). Irrespective of the value of p, the model presents first-order quantum 
transitions, when the long-range interaction is ferromagnetic, and second-order ones, when 



the long-range interaction is antiferromagnetic. Following Titvinidze and Japaridze 11], we 



have classified the intermediate phases, situated between saturated feromagnetic phases, as 

20 



the so-called quantum spin liquid phases, which are induced by the extended interaction. For 
p — 2, the global phase diagram has been obtained as a function of the interaction parameters 
and the critical surfaces and multicritical lines determined exactly. The critical exponents 
have been obtained and it has been verified that they satisfy the Rushbrook relation a+ 2(3+ 



7 = 2 (see e.g. ref. 28|) and the quantum hyperscaling relation 3Cj 2 — a = v (d + z) , and 



that the system presents a non-universal static and dynamic critical behavior. We have also 
shown that the free energy, close to first order transitions, satisfy the scaling form proposed 
by Continentino and Ferreira[l^]. 

Finally, we would like to point out that it is of great importance the existence of multiple 
first order quantum transitions and multicritical lines (triple and quadruple lines) , which 
we have shown exactly to exist in the model, since we believe that they are related to the 
different mechanisms from which the first order phase transitions are driven, as discusssed 



by Pfleiderer 



17j. 
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FIG. 1: Phase diagram for the quantum transitions as a function of the long-range interaction 
I/Ji for p = 2 and J2/J1 = 0.2,2.0. For J2/J1 = 0.2 there are three phases, one spin liquid phase 
(QSL-I) and two saturated ferromagnetic phases(SF), and for J2/J1 = 2.0, there are four phases, 
two spin liquid phases (QSL-I, QSL-II) and two saturated ferromagnetic phases(SF). The critical 
lines correspond to first-order phase transitions for I / J\ > and to second-order phase transitions 
for //Ji < 0. 
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I/J 1 

FIG. 2: Phase diagram for the quantum transitions as a function of the long-range interaction 
I/J\ for p = 2 and J2/J1 = —0.2, —2.0. For J2/J1 = —0.2 there are three phases, one spin liquid 
phase (QSL-I) and two saturated ferromagnetic phases(SF), and for J2/J1 = —2.0, there are four 
phases, two spin liquid phases (QSL-I, QSL-II) and two saturated ferromagnetic phases(SF). The 
critical lines correspond to first-order phase transitions for I / J\ > and to second-order phase 
transitions for I / J\ < 0. 
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FIG. 3: Magnetization M z as a function of h/J±, for p = 2, J2/J1 = 2.0 and different values I / J\. 
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FIG. 4: Isothermal susceptibility \t as a function of h/J±, for p = 2, J2/J1 = 2.0 and different 
values I/J\. 
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FIG. 5: Global phase diagram for p = 2 and J2/J1 > 0, as a function of h/J\ and I / J\. 
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FIG. 6: Global phase diagram for p = 2 and Jij J\ < 0, as a function of h/J\ and I / J\. 
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FIG. 7: Phase diagram for the quantum transitions as a function of the long-range interaction 
I/J± for p — ► oo and J2/J1 = 27/23, 2.0. For J2/J1 = 27/23 there are three phases, one spin liquid 
phase (QSL-I) and two saturated ferromagnetic phases(SF), and for J2/J1 = 2.0, there are four 
phases, two spin liquid phases (QSL-I, QSL-II) and two saturated ferromagnetic phases(SF). The 
critical lines correspond to first-order phase transitions for I / J\ > and to second-order phase 
transitions for I / J\ < 0. 
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FIG. 8: Phase diagram for the quantum transitions as a function of the long-range interaction 
I/Ji for p — ► oo and J2/J1 = —0.2, —1/11. For J2/J1 = —0.2 there are four phases, two spin liquid 
phases (QSL-I,QSL-II) and two saturated ferromagnetic phases(SF), and for J2/J1 = —1/11, there 
are three phases, one spin liquid phase (QSL) and two saturated ferromagnetic phases(SF). The 
critical lines correspond to first-order phase transitions for I / J\ > and to second-order phase 
transitions for I / J\ < 0. 
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FIG. 9: Phase diagram for the quantum transitions as a function of the long-range interaction I / J\ 
for p — ► oo and J^j ' J\ = —2.5, —1.7370... For J2/J1 = —2.5 there are four phases, two spin liquid 
phases (QSL-I,QSL-II) and two saturated ferromagnetic phases(SF), and for J2/J1 = —1.7370.., 
there are four phases, two spin liquid phases (QSL-I,QSL-II)) and two saturated ferromagnetic 
phases(SF) four phases which coexist in the quadruple point localized at h/J\ = 0, I / J\ = 1.7798... 
The functional of the free energy for this case is shown in Fig. [TT1 The critical lines correspond to 
first-order phase transitions for I / J\ > and to second-order phase transitions for I j J\ < 0. 



32 




FIG. 10: Phase diagram for the quantum transitions as a function of the long-range interaction 
I/J\ for p — > oo and J2/J1 = —3.5, —3.0. For J2/J1 = —3.5 there are four phases, two spin liquid 
phases (QSL-I,QSL-II) and two saturated ferromagnetic phases(SF), and for J2/J1 = —3.0, there 
are three phases, one spin liquid phase (QSL-I) and two saturated ferromagnetic phases(SF).The 
critical lines correspond to first-order phase transitions for I / J\ > and to second-order phase 
transitions for I / J\ < 0. 
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FIG. 11: Functional of the free energy / as a function of the magnetization M z for p — ► oo, 
h/Ji = 0, I/Ji = 1.7798... and J 2 /Ji = -1.7370... at the quadruple point. 
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FIG. 12: Phase diagram for the quantum transitions as a function of J2/J1 for p = 3 and // J\ = 0. 
All the transitions are of second-order and there are five phases, three spin liquid phase (QSL-I, 
QSL-II, QSL-III) and two saturated ferromagnetic phases (SF). 



35 




FIG. 13: Phase diagram for the quantum transitions as a function of J2/J1 for p = 4 and I / J\ = 0. 
All the transitions are of second-order and there are five phases, three spin liquid phase (QSL-I, 
QSL-II, QSL-IV) and two saturated ferromagnetic phases (SF). 



36 




FIG. 14: Static correlation function (SjSj +r )as a function of r, for p = 2, h/J± = 2.0, J2/J1 = 2.0 
and 1 1 J\ = —1.0. The inset shows the power law decay of the correlation characteristic of the 
QSL-I phase. 



37 




FIG. 15: Static correlation function (SjSj +r )as a function of r, for p = 2, h/J\ = —1.0, J2/J1 = 2.0 
and I/Ji = —1.0. The inset shows the power law decay of the correlation characteristic of QSL-II 
phase with the oscillatory modulation f(r). 
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FIG. 16: Static correlation function (SJ Sj +r )as a function of r, for p = 3, hj J\ = —0.15, J2/J1 = 
3.0 and I / J\ = 0. The inset shows the power law decay of the correlation characteristic of QSL-III 
phase with the oscillatory modulation f{r). 
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FIG. 17: Static correlation function (SjSj +r )as a function of r, for p = 4, h/ J\ = —0.6, J2/J1 = 1.6 
and //Ji = 0. The inset shows the power law decay of the correlation characteristic of QSL-IV 
phase with the oscillatory modulation f(r). 
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